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TWO-DIMENSIONAL  SYMMETRICAL  SUPERSONIC  NOZZLES 
By  J.  C.  Eward  and  Lawrence  R.  Marcus 


SUMMARY 

Auxiliary  boundary  conditions  are  derived  to  assure  continuity  of 
wall  curvature  in  applying  the  method  of  characteristics  to  the  design 
of  two-dimensional  symmetrical  supersonic  nozzles.  An  illustrative 
example  is  included. 


•  INTRODUCTION 

The  method  of  characteristics  was  first  applied  by  Prandtl  and 
Busemann  to  the  evaluation  of  supersonic  flows  in  1929  (reference  l) . 
Since  then^  a  number  of  expositions  (for  example^  see  references  2^  3^ 
and  4)  have  described  various  applications  of  characteristics  to  the 
design  of  two-dimensional  supersonic  nozzles.  The  problem  considered 
was  to  determine  wall  contours  that  would  transform  a  imiform  or  source 
flow  usually  at  Mach  number  of  unity  to  a  uniform  shock-free  flow  at 
some  higher  Mach  number. 

The  wall  contours  so  determined  usually  contain  discontinuities 
in  curvature  which  are  tolerable  for  fixed  geometry  designs.  In  varia¬ 
ble  Mach  number  flexible -walled  nozzles^  however^  a  discontinuity  in 
curvature  implies  a  multivalued  stress  at  a  pointy  which  in  turn  can 
only  occur  if  there  is  a  restraining  moment.  The  design  of  flexible- 
walled  supersonic  wind  tunnels  usually  does  not  include  provision  for 
such  a  restraining  moment  and  hence  the  wall  fairs  the  curvature  across 
the  discontinuity.  If  the  characteristic  diagram  for  the  tunnel 
requires  a  discontinuity^  then  imperfections  in  the  actual  flow  must 
appear . 

Even  if  a  rigid-walled  nozzle  is  used;  the  probability  that  the 
theoretical  gas  flow  will  occur  is  slight.  Changes  in  the  ratio  of 
specific  heats ;  or  any  condensation  of  the  gas  components;  or  imper¬ 
fections  in  the  corrections  for  boundary  layer  can  shift  the  character¬ 
istic  pattern  along  the  nozzle  wall.  The  seriousness  of  such  a  shift 
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Is  minimized  if  the  curvature  is  continuous  and  small.  Here  again  noz¬ 
zles  with  inherent  continuity  of  curvature  from  the  characteristic 
pattern  are  desired. 

A  plot  of  wall  angles  for  a  conventional  supersonic  nozzle  Including 
a  discontinuity  of  curvatixre  at  the  point  of  inflection  is  presented  in 
figure  1.  Also  shown  is  a  schematic  plot  on  which  the  discontinuity  in 
curvature  has  been  eliminated.  The  auxiliary  restrictions  which  must 
be  imposed  in  the  design  of  two-dimensional  supersonic  nozzles  to  elim¬ 
inate  discontinuity  of  wall  curvature  are  presented  herein.  The  work 
was  completed  during  May,  1951,  at  the  MCA  Lewis  laboratory. 


ANALYSIS  AND  DISCUSSION 

The  problem  is  to  determine  the  aTixiliary  conditions  which  must  be 
satisfied  to  eliminate  discontinuities  of  wall  curvature  in  the  design 
of  supersonic  nozzles.  In  conventional  designs,  a  discontinuity  occurs 
at  the  geometric  point  of  inflection  and  also  at  the  start  of  the  test 
s  ection . 

For  convenience,  two  sets  of  characteristics  ilr"^  and  T|f“  will  be 
utilized.  In  brief, 

t'*’  +  t”  =  ''  (l) 

is  the  Prandtl -Meyer  expansion  angle  tabulated  as  a  function  of  Mach 
number  in  reference  5.  The  quantity 

T(r+  -  \[f“  =  ^  (2) 

is  the  local  flow  deflection  angle  with  respect  to  the  axis.  Either 
or  \|f”  characteristics  represent  Mach  lines  in  the  flow.  The 
set  (fig.  2)  originates  on  the  contour  wall  and  moves  from  left  to 
right  as  the  characteristic  is  traced  from  the  wall  to  the  tunnel  center 
line.  The  f”  set  is  traced  from  left  to  ri^t  in  projecting  from  the 
axis  to  the  contour  streamline.  Thus ^  at  the  wall;  is  always  lar¬ 
ger  than  or  equal  to  . 

The  characteristic  point  is  that  place  on  the  wall  where  first 

assumes  the  test-section  value.  The  point  of  inflection  occurs  at  the 
position  of  maximum  wall  angle . 

The  curvature  along  a  streamline  may  be  expressed  mathematically 
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vhere  C  is  the  curvature  and  s  is  the  arc  length 
curvature  requires  that 


Continuity  of 


ds 


d\(f‘ 

ds 


d\|f’ 

ds 


(2a) 


be  continuous 


At  the  characteristic  point  and  downstream  therefrom^  is  con¬ 
stant  at  the  test-section  value.  The  curvature  is  then 


ds 


d\{r 

ds 


(3) 


downstream  of  the  characteristic  point.  The  point  of  inflection^  on 
the  other  hand^  has  zero  curvature  so  that  at  the  point  of  inflection 


d-ijf" 

ds 


d^ 

ds 


Thus^  if  the  characteristic  point  coincided  with  the  point  of  inflection^ 
continuity  of  curvature  would  require  that  the  local  values  of 

^  =  0.  A  zero  value  of  d\|/"/ds  at  the  point  of  inflection 
ds  ds  ds  ' 

results  only  if  di(r"‘'/ds  vanishes  at  some  point  on  the  wall  farther  upstream. 
Hence,  at  that  upstream  point,  d0/ds  =  -d-\|r“/ds.  Because  d'\|r“/ds  must  be 
positive  to  avoid  compression  waves,  d0/ds  would  then  be  negative  and  the 
nozzle  would  contain  at  least  two  inflection  points.  Rather  odd-shaped 
walls  therefore  result  which  are  undesirable  for  flexible  nozzles.  Thus 
the  characteristic  point  should  not  coincide  with  the  point  of  inflection 
if  the  wall  curvature  is  to  be  continuous,  (This  coincidence  does  occur 
on  most  conventional  nozzle  designs.)  Likewise,  by  similar  arguments,  the 
point  of  inflection  should  not  be  downstream  of  the  characteristic  point. 

The  remaining  possibility  is  for  the  point  of  inflection  to  be 
upstream  of  the  characteristic  point.  No  difficulty  is  experienced  in 
achieving  continuity  of  curvature  at  the  point  of  inflection  so  long  as 
diir“^  _  dilf” 

__  and  are  continuous  functions.  Continuity  can  be  achieved  at 

the  characteristic  point  also  if  approaches  zero  from  the  upstream 

side.  The  curvature  would  then  blend  smoothly  with  its  downstream  value 


given  by  equation  (3)  , 
to  its  final  value. 


Thus  t  is  tangent  at  the  characteristic  point 


Vanishing  of  at  the  characteristic  point  infers  a  vanishing 

at  the  start  of  the  test  section  (fig.  2).  Hence  the  discontinuity 

of  curvature  at  the  test -section  wall  junction  is  likewise  eliminated. 
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Thus  the  auxiliary  conditions  to  assure  continuity  of  wall  curva¬ 
ture  are: 

(1)  The  geometric  point  of  inflection  =  Oj  must  he  upstream 

of  the  characteristic  point. 

(2)  The  derivatives  and  must  he  continuous  functions, 

^  '  ds  ds 

and  as  a  conseq^uence  must  he  tangent  to  its  ultimate  test-section 

value  at  the  characteristic  point.  The  continuity  of  and 

ds  ds 

also  infers  continuity  of  Mach  number  gradient  along  the  axis . 

Nozzles  with  continuous  wall  curvature  will  always  he  longer  than  the 
shortest  possible  nozzle  for  a  given  Mach  number  and  test-section  size. 


An  exairiple  follows  for  which  the  wall  curvature  is  continuous. 

Plots  of  0,  and  Mach  number  distributions  for  an  M  =  4  nozzle 

are  presented  in  figure  3.  The  flow  was  expanded  from  a  uniform  Mach 
number  of  unity  at  the  throat  in  the  usual  manner  (reference  2)  to  the 
wall  point  of  inflection,  where  the  curves  of  and  \(f“  have  the 

same  slope.  The  \[f'^  curve  was  then  faired  from  tangency  at  the  point 
of  inflection  to  the  constant-test-section  value  in  a  horizontal  tangent 
at  the  characteristic  point.  The  remaining  values  of  i|f“  and  G  were 
determined  by  the  characteristic  plots  for  the  nozzle.  As  a  consequence 
of  the  fairing  at  the  characteristic  point,  0  and  ijr"  approach 

horizontal  tangents  at  the  test-section  junction  as  does  the  wall  Mach 
number . 

In  practice,  nozzles  of  predetermined  length  are  often  desired. 
Precise  estimation  of  the  nozzle  length  for  a  given  test-section  size 
is  difficult  when  only  the  characteristic  point  and  final  Mach  number 
are  specified.  An  alternative  design  procedure  is  therefore  useful. 

The  flow  in  the  vicinity  of  the  nozzle  throat  can  be  obtained  by  charac¬ 
teristics  in  the  usual  manner.  A  curve  of  either  or  Mach  n*umber  is 

thus  obtained  along  the  axis  which  properly  matches  the  initial  super¬ 
sonic  flow  to  the  sonic  flow  at  the  throat.  From  the  design  Mach  number 
and  test-section  size,  the  upstream  vertex  of  the  test  diamond  is  then 
located.  Continuity  of  dijr+Zds  and  dijr^/ds  along  the  nozzle  contour 
implies  continuity  of  dafr'^/ds  and  also  dM/dx  along  the  axis.  Either 
the  value  or  the  Mach  number  distribution  along  the  axis  should  then 

be  faired  in  a  smooth  curve  to  a  horizontal  tangent  at  the  upstream  vertex 
of  the  test  diamond.  Working  from  the  axis  the  remaining  wall  contours 
are  determined  by  the  characteristic  plots  to  give  a  nozzle  of  the 
required  length  and  size. 


CONCLUSIONS 

Auxiliary  boundary  conditions  were  derived  to  assure  continuity  of 
wall  curvature  in  the  application  of  the  method  of  characteristics  to 
two-dimensional  symmetrical  supersonic  nozzles.  These  conditions  are: 
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1.  The  point  of  inflection  must  be  upstream  of  the  chaxact eristic 
point . 


2.  The  first  derivatives  of  the  characteristic  values 


^  and 
ds 


must  be  continuous  functions .  Hence  must  be  tangent  to  its 

test-section  value  at  the  characteristic  point.  This  latter  restriction 
infers  that  the  wall  slope  approaches  a  minimum  or  horizontal  tangent 
(that  is^^zero  curvature)  at  the  test-section  junction.  Also  the  axial 
Mach  number  gradient  must  be  continuous  and  approach  zero  at  the  upstream 
vertex  of  the  test  diamond. 


Lewis  Flight  Propulsion  Laboratory 
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Cleveland^  Ohio^  October  16^  1951 
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Axial  distance^  x 


Figure  3,  -  Variation  of  flow  parameters  for  a  Mach  numher  4  tvo- dimensional 

supersonic  nozzle. 
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